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The fermion determinant in an instanton background for a quark field of arbitrary mass is de- 
termined exactly using an efficient numerical method to evaluate the determinant of a partial wave 
radial differential operator. The bare sum over partial waves is divergent but can be renormalized in 
the minimal subtraction scheme using the result of WKB analysis of the large partial wave contri- 
bution. Previously, only a few leading terms in the extreme small and large mass limits were known 
for the corresponding effective action. Our approach works for any quark mass and interpolates 
smoothly between the analytically known small and large mass expansions. 



To study instanton-related physics in QCD it is of 
fundamental importance to determine the one-loop tun- 
nelling amplitude given by the Euclidean one-loop effec- 
tive action in an instanton background Q] . This quantity 
involves the fermion determinant in a nontrivial gauge 
configuration and so is also of interest for unqucnching 
effects in dynamical quark simulations in lattice QCD 2] . 
In the 1970s, 't Hooft 0] calculated analytically this one- 
loop amplitude for massless scalar or quark fields, but 
this exact calculation is no longer possible if the fields 
have nonzero mass. The small mass limit was extended 
further by Carlitz and Creamer 4 : ], and by Kwon et al. 
0. This small mass limit is also closely related to the 
study of zero modes and spectra of the Dirac operator, 
which have been investigated extensively recently in lat- 
tice QCD 0. In the other extreme, the large mass limit 
is naturally obtained from the Schwinger-DeWitt (or the 
heat-kernel) expansion pllflltlUi] within the proper-time 
representation of the effective action. For phenomeno- 
logical applications 0, and also for the extrapolation of 
lattice results 0| obtained at unphysically large quark 
masses to lower physical masses, it is important to be able 
to connect the large and small mass regimes. In this pa- 
per, we provide the first precise bridge between these two 
extremes of small and large mass, by presenting a method 
which computes the exact effective action in an instanton 
background for any value of the quark mass. The result- 
ing mass dependence interpolates smoothly between the 
explicitly known results at the opposite ends. 

Our computational procedure is based on an efficient 
way to compute one-dimensional determinants, combined 
with the fact that an instanton background is sufficiently 
symmetric that the corresponding effective action can be 
reduced to a sum over partial waves (see details below). 
This bare sum over partial waves is of course divergent, 
and the nontrivial part of the calculation is to renor- 



malize this divergent sum in a physically unambiguous 
manner. We use the minimal subtraction scheme 0,0 
and a WKB analysis > which gives the exact countert- 
erms, including finite parts. The result is a simple finite 
expression for the renormalized one-loop effective action 
[see (01 below] which can be evaluated numerically. 

Due to a hidden supersymmetry 0, 0], the spinor 
Dirac operator in an instanton background has the same 
spectrum (except for zero modes and an overall multi- 
plicity factor of 4) as that of the scalar Klein-Gordon 
operator in the same instanton background. This im- 
plies that the one-loop effective action of a Dirac spinor 
field of mass m (and isospin i), Tf en (A;m), is directly 
related to the corresponding scalar effective action (fo r 
a complex scalar of mass m and isospin i ) by 0,0, Il4| : 

rf cn (A;m) = -|ln(=£) - 2Tf cn (A;m), where /x is a 
renormalization scale. The logarithmic contribution cor- 
responds to the existence of a zero eigenvalue in the spec- 
trum of the Dirac operator for a single instanton back- 
ground. This relationship (which is special to a self-dual 
background) has the important consequence that it is suf- 
ficient to consider the scalar effective action Tf cn (A;m) 
to learn also about the corresponding fermion effective 
action T^(A; m), for any mass m. 

As in we consider an SU(2) single instanton of scale 
p : A^(x) — v ^ 2 a + p 2 V ■ Rather than the scalar effective 

action Tf en (A;m), it is convenient to consider rf cn (mp), 
a function of mp only, defined by 



Tf ca (A; m) = - ln(pp) + ff cn (mp). 



(1) 



The factor g in (JIJ is determined by the one-loop (3 func- 
tion. We can then concentrate on the mp dependence of 
ff cn (mp), and we can set the instanton scale p = 1. In 
the small and large mass limits, it is known that rf cn (m) 
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behaves as 0, 
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where a(|) ~ 0.145873, and 7 ~ 0.5772... is Euler's 
constant. These extreme small and large mass limits 
are shown as dashed curves in Figure ^ There is clearly 
a significant gap preventing extrapolation between the 
small and large mass limits, and for mass between 0.5 
and 1 (in units of 1/p) neither of these approximations is 
particularly accurate. In this paper we provide a simple 
numerical procedure to determine Tf cn (m), and hence the 
fcrmion determinant, for any value of the mass to, not 
just asymptotically small or large masses. Our result is 
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where L is a large integer. The sum over partial- wave 
contributions Tf is done numerically, as described below, 
in Eq. (|13|l . The other terms in are renormalization 
terms, computed using minimal subtraction and WKB. 
The renormalized effective action Tf on (m) in © is finite, 
converges for large L, and can be computed for any mass 
771. Figure 1 shows that our numerical results provide a 
very precise interpolation between the extreme small and 
large mass limits in (j2J|. 
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FIG. 1: Plot of our numercial results for ff en (m) from 10, 
compared with the analytic extreme small and large mass 
limits [dashed curves] from 0. The dots are numerical data 
points from ©, and the solid line is a fit through these points. 



To derive © , note that the regularized one-loop scalar 
effective action has the proper-time representation 



r x (T„M = - / ^(e- 



< J d'j- ir(r < 
ds 
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From this one obtains the renormalized effective action, 
in the minimal subtraction scheme 0, H : 



r? en (A;m) 



lim 



r A (A;TO)--t-ln( ^ 
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Recall that the massive Klein-Gordon (KG) operator, 
— D 2 + to 2 , for scalars of isospin i in the instanton back- 
ground (with scale p = 1) reduces to the radial form Q 



M {hj) = ----- 



a 2 

dr 2 
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where Z = 0, 



1, 



. . ; j = and with degeneracy 

, (3) factor dnf\ = (2Z + l)(2j + 1). Without the instanton 



background, the free Klein-Gordon operator is 
M?r = -—-~- + ^+V+rn 2 



(7) 



with degeneracy factor = (2Z+1) 2 . Thus, it is natural 
[l2| to combine the partial waves (l,j = I + | ) and (/ + 
j = /), which have common degeneracy factor (21 + 
1) (2Z + 2), so that the sum over I and j reduces to a 
sum over I. This sum is, however, divergent and we must 
define a consistent regularization and renormalization. 
To this end, we use the proper-time regularization in Q , 
and split the I sum as follows: 

L 00 

r s A (A;m) = ^ rf jW (^;m) + ^rf >w (^;m)(8) 
;=o,§,... 

where L is a large but finite integer. In the first sum, 
which is finite, the cutoff A is irrelevant and we can 
use a numerical method [described below] to evaluate 
the sum. For the second sum, we use a combination of 
the WKB approximation, which is good for large I, and 
Euler-Maclaurin summation to perform the sum. 
To begin, we consider the first sum in ©: 
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These one-dimensional determinants can be computed ef- 
ficiently using the following result El El El El El : 
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Suppose Mi and M 2 are two second order ordinary dif- 
ferential operators on the interval r £ [0, 00), with Dirich- 
let boundary conditions. Then the ratio of the determi- 
nants of A4i and M. 2 is given by 
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where ipi( r ) = li 2) satisfies the initial value problem 

Mi^i(r) = ;V)i(r = 0) = 0;V-(^ = 0) = l (11) 

Since an initial value problem is very simple to solve nu- 
merically, this theorem provides an efficient way to com- 
pute the determinant of an ordinary differential operator. 
Note in particular that no direct information about the 
spectrum (either bound or continuum states, or phase 
shifts) is required in order to compute the determinant. 

We can simplify the numerical computation further be- 
cause for the free KG operator, M f ^ c , the solution to ltTT|) 



L (mr) 



is just the modified Bessel function: V , m a ( r ) = — 
Since this grows exponentially fast at large r, this should 
also be true of the numerical solutions of Hll|) for the 
instanton operators, -Mnj), in ©. Thus, it is numeri- 
cally better to consider the ODE satisfied by the ratio (a 
similar idea was used by Baacke et al in their analysis of 
metastable vacuum decay [2jj) : T( Zj )(r) = ^tll^ ■ In 

fact, since we are ultimately interested in the logarithm of 
the determinant, it is more convenient (and more numer- 
ically stable) to consider S^j) = hiTnj\, which satisfies 



d 2 S, 



{1,0) 



1 



dr 2 \ dr 

_ 4(j-Z)(j + Z + l) 
r 2 + 1 



2m 



I! 2i+i( mr )\ dS (i,i) 



hi+i(mr) J dr 
3 



(r 2 



1 



(12) 



where Snj-\(r) satisfies the initial conditions Snj\(0) = 
and S'n jJO) = 0. Thus, the first, finite, sum in © for 
the bare effective action can be evaluated as 



Tf(A;m) = J2 (2? + 1)(2/ + 2) 



1=0, i, 



tO,i 



x |5 (M+ i)(r = 00) + S (l+ i tl) (r = 00)} .(13) 

The large r values of S^^ + i^(r) and S^ l+ i «(r) can be 
extracted with excellent precision. In fact, the asymp- 
totic values of ; + i)(r) and S^ + i ^(r) very nearly can- 
cel one another, and for a given mass, as a function of Z, 
the combination [S^ ^i) (r = 00) + S^ + i ^(r — 00)] falls 
off in magnitude like j - see Figure |2 Thus, the sum in 
(|13J) will have terms going like, L 2 , L and InL, as well as 
terms finite and vanishing for large L. We now show that 
these potentially divergent terms are exactly canceled by 
terms in the second sum in 
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FIG. 2: Plot of the I dependence of P(l) = S (l l+ iJr = 00) + 
q (r = 00), for m = 1. The sum of these two asymptotic 
values is O(j) for large I. 



The second sum in (ISJ can be analyzed using the Euler- 
Maclaurin method [2l| as follows. Write 



x poo r I . 

J2Tl (l) (A;m)= ds --( 
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Using WKB, which is good for large Z, we can write 01 
/(*)= / <lr ! >, .//(-/•) i CL- r .) 



/ dr £ / ; ( S ,r) 



where for each Z, fi(s,r) has a local expansion in terms 
of the Langer-modified potential j2^ 
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■(16) 



The first three orders of the WKB approximation for /; 
were computed in \j/2}. We find fi(s,r) = {21 + 1)(2Z + 
2)[/ (M+ i ) (s,r) + / (i+ i :i) (s,r)], where to first order 



and to second order 
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The sum over Z in (|15|l can now be performed using the 
Euler-Maclaurin expansion [2l| : 

00 />oo -1 -1 

£ /i =2/ ^ /(0 - o /(L) " 24 /,(i) + • ' ■ (19) 

For each order of the WKB expansion, the I dependence 
is such that the Z integral in l|19|) can be done exactly. 



4 



Then in the large L limit, the r and s integrals in 1|15|) 
and l|14fl can also be done, leading to |2^| 



r A,(0 



l=L-i 



^ lnA 

6 



2L 2 + AL 



+ 



127 
~72 



3 



■ In 2 



m 
~~2 



m 2 In 2 - 



m 
~2~ 



1 
6 

In to 



(20) 



Only the first two orders of the WKB expansion are re- 
quired for this result; all higher orders contribute terms 
vanishing in the large L limit. It is important to iden- 
tify the physical role of the various terms in l|2(J|l . The 
first term is the expected | In A term canceled in (JSJ . 
The next three terms give the quadratic, linear and loga- 
rithmic divergences in L which cancel the corresponding 
divergences in the first sum in JSJ, which were found in 
our numerical data - see Fig. [2] It is a highly nontrivial 
check on the WKB computation that these terms have 
the correct coefficients to cancel these divergences. Com- 
bining (|2"U|> and (|T3|) . we obtain our result ©. 

As an interesting application, our result © provides 
a very simple analytic computation of 't Hooft's leading 
small mass result. When to = 0, the numerical integra- 
tion can be done exactly and one finds that |23| 



S (l,l+k)( r = °°) + S Q+i,l)( r = °°) 
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Then it follows from @ that 
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- = 0.145873. 
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